Condensate fraction in neutron matter 
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We study the Bose-Einstein condensation of fermionic pairs in the uniform neutron matter by 
using the concept of the off-diagonal long-range order of the two-body density matrix of the system. 
We derive explicit formulas for the condensate density peon and the condensate fraction peon I p 
in terms of the scaled pairing energy gap A/ef where is the Fermi energy. We calculate the 
condensate fraction peon I p as a function of the density p by using previously obtained results for 
the pairing gap A. We find the maximum condensate fraction {peon/ p)max = 0.42 at the density 
p = 5.3 • 10~* fm~^, which corresponds to the Fermi wave number kp = 0.25 fm~^. 
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PACS numbers: 21.65.Cd, 03.75.Ss 

The concept of off-diagonal long-range order (ODLRO) 
was introduced by Penrose and Onsager [1| to determine 
the condensate fraction of a generic bosonic many-body 
system. Yang extended this powerful idea to the case 
of fermionic many-body systems. It is now established 
that at zero temperature the condensate fraction of liq- 
uid ^He is below 10% Q, while for dilute and ultracold 
bosonic alkali-metal atoms it can reach 100% 0]. Few 
years ago the condensate fraction has been calculated for 
fermionic atoms [5|, |6j in the crossover from the Bardeen- 
Cooper-Schrieffer (BCS) state of Cooper Fermi pairs to 
the Bose-Einstein condensate (BEC) of molecular dimers 
at zero temperature. These theoretical predictions are in 
quite good agreement with the data obtained in two ex- 
periments with ultra-cold two hyperfine component 
Fermi vapors of ^Li atoms. Very recently, the conden- 
sate fraction of dilute neutron matter, which is predicted 
to fill the crust of neutron stars [31 , has been calculated 
by Wlazlowski and Magierski [lO, [lH. Unfortunately, 
these finite-temperature Path Integral Monte Carlo sim- 
ulations are performed only at p = 0.003 fm"'^, and their 
trend suggests a condensate fraction of about 35% at zero 
temperature [Tll |. 

In this paper we present a zero-temperature system- 
atic study of the condensate fraction of neutron matter 
as a function of the density p. First we derive explicit 
analytical formulas for the condensate density pcon and 
the condensate fraction Pconj P as a function of the pair- 
ing energy gap A in units of the Fermi energy tp- With 
the help of these expressions we calculate the condensate 
fraction PconI P vs p by using previously obtained results 
for the paring gap A [ij, [ij] . 

The shifted Hamiltonian of the uniform neutron matter 
is given by 

H' = [d'r J2 rpUr)(-^y'-Ai'Ar) (1) 
+ y d\ d'r' i^+{r) ^+{r') V{r - r') ^;(r') V't(r) , 

where i^cri^) is the field operator that annihilates a neu- 
tron of spin a in the position r, while ipai^) creates a 



neutron of spin a in r. Here V{r — r') is the nucleon- 
nucleon potential in the ^Sq channel. The ground-state 
average of the number of neutrons reads 

j d^v J2 (^^W Mr))= I rf'r (p(r)) . (2) 

This total number A'' is fixed by the chemical potential 
fi which appears in Eq. ([T]). 

As shown by Yang 0], for a Fermi system ODLRO 
means that the two-body density matrix factorizes as fol- 
lows 
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in the limit wherein both unprimed coordinates approach 
an infinite distance from the primed coordinates. The 
largest eigenvalue No of the two-body density matrix 
gives the number of neutron pairs which have their cen- 
ter of mass with zero linear momentum. This condensed 
number of pairs is given by 5, 6] 



Ar,„„ = 2 / d^n d^V2 |( V;(ri) ^t(r2))r • (4) 



Within the Hartree-Fock-Bogoliubov approach [ij], 
the shifted Hamiltonian density Q can be diagonalized 
by using the Bogoliubov-Valatin representation of the 
field operator tpa- (r) in terms of the anticommuting quasi- 
particle Bogoliubov operators 6ko- with amplitudes Uk 
and and quasi-particle energy i^k- After minimization 
of the resulting quadratic Hamiltonian one finds familiar 
expressions for these quantities (l3 |: 
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FIG. 1: Condensate fraction peon/ p as a function of the scaled 
pairing gap A/ef- The sohd hne is obtained by using Eq. 
(|14p . The dashed hne is obtained by using Eq. (|15|l . valid in 
the regime A/ef ^ 1- 



where ek is the single-particle energy and Ak is the pair- 
ing gap, which satisfies the integral equation 



A, 



(7) 



In the Hartree-Fock-Bogoliubov approach [ij] the single- 
particle energy is given by ek = ek + U\^, where 
ek = h^k"^ / {2m) is the bare single-particle energy and 
[/k = X^qC^i q) is the Hartree-Fock potential. The 

equation for the average number of neutrons is easily ob- 
tained from Eq. ([2]) as 
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(8) 



Finally, from Eq. (|4]) one finds that the condensate num- 
ber of paired neutrons is given by 
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The numerical solution of Eqs. ([7]), ([S]), and ^ is 
a hard computational task. However, the problem can 
be simplified under some reliable assumptions In 
Eqs. ^ we adopt the approximations ek — ek -I- Ukp 
and /i ~ e^F + Ukp, where kp — (Stt^p)^/'^ is the Fermi 
wavenumber of the non-interacting Fermi gas with den- 
sity p. In this way we get Ck — A* — h^k"^ /{2m) — ep where 
eF = ^kp = fi.^(37r^p)^/'^/(2m) is the Fermi energy of the 
non-interacting Fermi gas. Notice that this is formally 
equivalent of setting directly Ck — Ck and ~ e^? in Eqs. 
©. In addition, we set Ak — A^^ = A. Under these as- 
sumptions, in the three-dimensional continuum limit the 
gap equation ([7]) of the neutron matter becomes 
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Moreover, from the number equation ([8]) we find the total 
density as 
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where I{x) is the monotonic function 
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In a similar way from Eq. (|9|) we get the condensate 
density of the neutron-neutron pair 
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The condensate fraction follows immediately as the ratio 
between pcon , given by Eq. (fT3|) and p, given by Eq. (|lip , 
namely 



A ^ 
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Notice that this fraction is expressed only in terms of the 
ratio ep/A. In the deep BCS regime where A/cp ^ 

V2 {ep/A)^^'^ and 



1 then Jep/A+ ^1 + 4JA2 



I {ep/A) {2/3){ep/A) 
densate fraction becomes 
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and consequently the con- 
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In Fig. [T] we plot the condensate fraction Pcon/p as a 
function of the scaled pairing gap A/ep. The solid line 
is obtained by using Eq. (|14p . while the dashed line is 
obtained by using Eq. (|15p . The figure shows that for 
A/e_F < 0.2 the two curves practically coincide. 

To estimate the condensate fraction Pcon/ P in neutron 
matter on the basis of Eq. (fT^ it is necessary to ob- 
tain A/ep from Eq. (jlOp or, if possible, directly from 
Eq. ([7]). This problem has been investigated by various 
authors [H HI |il-il| . For very dilute neutron matter, 
the nucleon-nucleon potential V{y—y') in the ^Sq channel 
can be parametrized in terms of only the s-wave scatter- 
ing length tts = —18.5 fm [ij]- Indeed, in the regime 
/cfIosI ^ Ij the pairing gap is exponentially small and 
from Eq. (jTU]) one finds 
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It follows immediately from Eq. (|15p that in this regime 
of very- low density (BCS regime) the condensate fraction 
reads 
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FIG. 2: (Color online). Scaled pairing gap A/ep of neutron 
matter as a function of the Fermi wave number kp- The 
dashed line is obtained with Eq. (|16p : the filled squares with 
the dot-dashed line are the results obtained with the GSRS 
nuclear potential [13]; the filled circles with the solid line are 
the results obtained with the Argone V18 nuclear potential 

El. 



At higher densities the pairing gap A does not follow an 
exponential grow but instead it reaches a maximum and 
eventually goes to zero [H [l3, [H-ill- In Fig. [2 we 
report the scaled pairing gap A/ep as a function of the 
Fermi wave number kp. The solid line, obtained with Eq. 
(ITS)) , shows the exponential grow which is reliable only in 
the very low density BCS regime. The filled squares are 
the results obtained by Matsuo 112)] by solving Eq. ([7]) 
with the G3RS nuclear potential [23]; the filled circles 
are the results obtained by Gezerliz and Carlson fl^ with 
the Argone V18 nuclear potential [23]. The figure shows 
that the two nucleon-nucleon potentials give very similar 
results. 

Instead of solving Eq. pU)) we use the data of A/ep 
shown in Fig. [5] to derive the condensate fraction Pcon/p 
by using our Eq. 1^^. It is then physically relevant to 
plot the obtained condensate fraction Pcon/p as a func- 
tion of the density p = fc|,/(37r^). 

The results are shown in Fig. [3l where the horizontal 
axis is in logarithmic scale. At very low density p the 
neutron matter behaves like a quasi-ideal Fermi gas with 
weakly correlated Cooper pairs and the condensate frac- 
tion Pcon/p is exponentially small (BCS regime, dashed 
line). By increasing the density p the attractive tail of 
the neutron-neutron potential becomes relevant and the 
condensate fraction pcon/ P grows significantly. On the 
basis of the data obtained with the Argone ^18 potential 
(filled circles with dot-dashed line) , the maximum of the 
condensate fraction is {pcon/ p)max = 0.42 at the density 
p = 5.3 • 10~^ fm~'^ which corresponds to the Fermi wave 



number kp = 0.25 fm~^. By further increasing the den- 
sity p the repulsive core of the neutron-neutron potential 
plays an important role in destroying the correlation of 
Cooper pairs and the condensate fraction Pcon/ P slowly 



0.5 RTT 

0.45 - 
0.4- 
0.35 - 
Q. 0.3- 
o 0.25 - 
Q- 0.2- 
0.15 - 
0.1 - 
0.05 - 




10 10 10" 10 10"^ 



10" 10 



FIG. 3: (Color online). Condensate fraction pcon /p of neutron 
matter as a function of the density p. It is derived by using 
the data of A/ef of Fig. [2] and Eq. (fTi)) . 



goes to zero (filled squares with solid line). Remarkably, 
the results of Fig. [T] are fully consistent with the Monte 
Carlo value Pcon/ P — 0.35 at p = 0.003 fm^'^ one ex- 
tracts from the finite-temperature Path Integral Monte 
Carlo data of Wlazlowski and Magicrski 10, 11]. 

As discussed in the introduction, the properties of low- 
density neutron matter are important for the understand- 
ing of neutron star crusts 9] and also the exterior of 
neutron-rich nuclei [2^ . Our present investigation gives 
a contribution to this fascinating field of research: the 
condensate fraction of neutron matter reaches its maxi- 
mum value of about 40% in correspondence to the max- 
imum of the pairing gap at the density of about 5 • 10~^ 
fm~'^, which is much smaller than the nuclear satura- 
tion density 0.16 fm~^. It is important to stress that 
neutron-neutron pairs do not form a true bound state 
and the system cannot reach the unitarity limit (where 
there is the formation of the bound state) nor the deep 
BEC regime characterized by a gas of di-neutrons with 
100% condensate fraction [T^. Nevertheless, the maxi- 
mum value of the condensate fraction we have found for 
neutron matter is not too far from to the value of about 
70% calculated for the unitarity Fermi gas [5] . Indeed our 
results clearly show a BCS-quasiunitary-BCS crossover 
by increasing the neutron density. 

The author thanks Prof. Masayuki Matsuo for making 
available his numerical data and for useful suggestions. 
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